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Using density functional theory we investigate the evolution of the magnetic ground state of 
NbFe2 due to doping by Nb-excess and Fe-excess. We find that non-rigid-band effects, due to the 
contribution of Fe-d states to the density of states at the Fermi level are crucial to the evolution of 
the magnetic phase diagram. Furthermore, the influence of disorder is important to the development 
of ferromagnetism upon Nb doping. These findings give a framework in which to understand the 
evolution of the magnetic ground state in the temperature-doping phase diagram. We investigate the 
magnetic instabilities in NbFe2- We find that explicit calculation of the Lindhard function, Xo(q), 
indicates that the primary instability is to finite q antiferromagnetism driven by Fermi surface 
nesting. Total energy calculations indicate that q = antiferromagnetism is the ground state. We 
discuss the influence of competing q = and finite q instabilities on the presence of the non-Fermi 
liquid behavior in this material. 

PACS numbers: 75.10.Lp, 71.20.Be, 71.20.-b 



I. INTRODUCTION 

Materials that exhibit magnetic phase transitions that 
may be tuned to zero temperature are strong candidates 
for the formation of novel ordered states around the 
quantum critical point. These novel states rely upon the 
breakdown of assumptions that underpin the Fermi liquid 
theory of metals when the magnetic interactions become 
extended in space and time. Such a situation is possible 
at a magnetic quantum critical point and may lead to the 
formation of new phases of matter that are potentially 
technologically valuable. The formation of superconduc- 
tivity in cuprates and iron-pnictides near magnetic phase 
transitions are among the phenomena that motivate us 
to understand materials where magnetic interactions are 
potentially of crucial importance. 

The metallic compound NbFe2 has been of interest to 
researchers for decades due to its rich magnetic phase di- 
agram. Initially, stoichiometric NbFe2 was thought to be 
either paramagnetic or ferromagnetic^. However, sub- 
sequent NMR 3 and magnetization measurements^ have 
been interpreted as showing spin density wave (SDW) 
order with a Neel temperature T/v ~ 10K. 

Studies of the doping evolution of the phase diagram 
have shown that this spin density wave state is very sen- 
sitive to changes in stoichiometry. Studies of doping in 
Nbi_ y Fe2+ 2/ with either Nb-excess (y < 0) or Fe-excess 
(y > 0) have shown that stoichiometric NbFe2 becomes 
ferromagnetic upon either electron or hole doping 5 . Be- 
cause of the close proximity of the reported SDW and 
ferromagnetism, NbFe2 has long been considered a can- 
didate material for the coexistence of SDW and ferro- 
magnetic fluctuations^. 

Recent experimental work has further probed the 
behavior of the temperature-doping phase diagram of 
high purity NbFe2 via measurements of the tempera- 
ture dependence of resistivity, specific heat and magnetic 
susceptibility^. The resulting phase diagram is shown 
in Fig. [T] Importantly these refined measurements in- 
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FIG. 1: (Color online) Phase diagram for Nbi_ y Fe2+ y 7 . 
By adjusting the precise composition within a narrow range, 
NbFe2 can be tuned from ferromagnetism (FM) (y > 0.01) via 
an intermediate spin-density-wave (SDW) modulated state to 
a quantum critical point (QCP) (y ~ —0.015). 



dicate the presence of a quantum critical point as the 
reported SDW collapses when NbFe2 is tuned away from 
stoichiometry by Nb-excess. At this critical point, in- 
triguing non-Fermi liquid behavior in the material has 
been demonstrated. Wada et a/P had previously re- 
ported an enhanced electronic specific heat in the late 
1980's. This behav ior has been confirmed by recent 
measurement JUSE] that describe the non-Fermi liquid 
behavior in more detail. They show at low temperatures 
near the quantum critical point the presence of a log- 
arithmic divergence of the specific heat and a p ex T2 
dependence in the resistivity. 

Further NMR and magnetization studies have 
confirmed the presence of antiferromagnetism at 
stoichiometry^E^n^, Dl rt neutron diffraction^^ studies 
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have been unable to determine the magnetic order. 
Therefore, studying the electronic and magnetic prop- 
erties of this material with electronic structure calcula- 
tions may aid our understanding and has been addressed 
previousljPEl. Xotal 

energy calculations with the frozen 
core approximation^ suggested a paramagnetic ground 
state, while work with the LMTO method with shape 
approximations 16 suggested an antiferromagnetic ground 
state nearly degenerate with a paramagnetic state. 

Recent work by Subedi et al. has discussed in de- 
tail the magnetic interactions of local moment configu- 
rations within the unit celP^ of stoichiometric NbFe2- 
NbFe2 exhibits significant exchange splitting of the 3s 
core levels 18 . This indicates the presence of local 
moment fluctuations as is also found in the Fe-based 
superconductors^. Subedi et al. found that the mag- 
netism has itinerant character due to the dependence 
of the size of the moments on the magnetic configura- 
tion. Also, a ferrimagnetic ground state was found to 
be most favourable. Furthermore large moments ~ \[Ib 
were obtained, which far exceed the ~ 0.02/i^ inferred 
from experiment. In this study we further address the 
magnetic ground state of this material using the full po- 
tential methods in WIEN2k 20 . We consider the doping 
dependence of the ground state and discuss the magnetic 
interactions that may be key to the critical behavior in 
this material. Our findings for the ground state are used 
to discuss the driving forces behind the rich magnetic 
phase diagram and the role of magnetism in the intrigu- 
ing non-Fermi liquid behavior of this material. 

We begin our analysis by considering the electronic 
structure of stoichiometric NbFe2- All calculations are 
performed using the GGA-PBE 21 correlation functional 
of Perdew, Burke and Ernzerhof. We have used the 
experimental lattice parameters^. For stoichiometric 
NbFe 2 a = 4.8401(2) A and c=1.8963(6) A. All results 
presented were obtained using RK max = 8, G max = 16 
and IFFT factor=4.0. Here, RK max determines the ma- 
trix size (K max is the plane wave cut-off, R is the small- 
est of the atomic sphere radii), G max is the maximum 
wavevector in the charge density Fourier expansion and 
IFFT (Indices of the Fast Fourier Transfrom grid) de- 
termines the size of the mesh for the calculation of the 
exchange-correlation potential in the interstitial region. 
The Fermi surface and charge densities were calculated 
from a 43 x 43 x 23 grid of k-points. The radii of the muffin 
tins were 2.12ao for Fe and 2.15ao for Nb. 



II. CRYSTAL AND BONDING STRUCTURE 

NbFe2 crystallizes in the C14 Laves phase which has 
the MgZn2 hexagonal structure. The crystal structure is 
shown in Fig. [2j The Fe atoms form a layered structure 
of a kagome lattice (6/i sites) separated by Fe (2a sites) 
atoms centred on the line between alternate kagome tri- 
angles. The Nb atoms occupy the interstices in this Fe 
structure and lie slightly out of plane with respect to the 




FIG. 2: (Color online) Crystal Structure of NbFe 2 . Fe(6/i) 
red, Fe(2a) blue and Nb grey. 



Fe 2 a sites. The site symmetries have internal degrees 
of freedom: Nb at 4/(l/3/,2/3,x) and Fe at 2a(0,0,0) 
and 6/i(y,2y,3/4). In our calculation of stoichiometric 
NbFe2 we have relaxed these internal parameters result- 
ing in x=0.0657 and y= 0.1705. 

It is worth considering this crystal structure in further 
detail since its properties affect the bonding, doping im- 
purity location and the magnetic coupling. In Fig. [3ja) 
we show an extended view of the crystal structure along 
the c-axis. We can see that the blue Fe(2a) atoms only 
coordinate every second triangle of the kagome structure 
formed by the Fe(6/i) sites. To emphasise this we have 
separated the colors of the upper and lower kagome layers 
into red and green. The presence of the Fe(2a) atoms co- 
ordinating every second kagome triangle means that the 
bonding between the pairs of Fe(6/i) atoms in the kagome 
layers will not be the same. In Fig. |3|b) we show charge 
density contours within the kagome layer. As we can 
see the bonding between Fe(6/i) sites varies significantly. 
Effectively the presence of the Fe(2a) sites draws away 
charge density into the out of plane Fe(2a)-Fe(6/i) bonds, 
and starves the bonds in the kagome triangles that are 
coordinated with an Fe(2a) site. 

In Table [J we show the interatomic bond distances. 
Significantly, amongst the Fe sites, the nearest neigh- 
bour distances from Fe(2a) and Fe(6/i) sites are similar. 
Therefore upon Nb doping, the preferred dopant site is 
likely to be determined by the bonding network of the 
Fe cage rather than the volume available at the site. We 
investigate this via total energy calculations in the fol- 
lowing section. 

The density of states for non-magnetic NbFe2 is 
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FIG. 3: (Color online) Crystal Structure of NbFe 2 at left. 
View shown along the c-axis. The upper and lower kagome 
layers formed by the Fe(6/i) sites have been separated into red 
and green respectively. Fe(6/i) red (upper kagome), Fe(6/i) 
green (lower kagome), Fe(2a) blue and Nb grey. In the right 
figure we show charge density contours in the upper kagome 
plane. The charge density for kagome triangles coordinated 
with Fe(2a) sites is inequivalent to those with no coordination. 
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FIG. 4: (Color online) Non-magnetic calculated DOS for 
NbFe 2 . The Fermi level is set to eV. 



TABLE I: Interatomic bond distances for stoichiometric 
NbFe2- In each case we show the shortest distance between 
the two given sites. 



Sites 



Distance (A) 



Fe(2a)-Fe(6/i) 
Fe(6/i)-Fe(6/i) 
Fe(2a)-Fe(2a) 
Nb-Fe(2a) 
Nb-Fe(6/i) 
Nb-Nb 



2.42 
2.37 
3.95 
2.84 
2.81 
2.89 



shown in F ig. [4] and agrees with the results of previous 
studies^E^l We also show the partial density of states 
by atomic site. The partial density of states for Fe(2a)-d 
and Fe(6h)-d sites possess similar structure, which is con- 
sistent with their strong bonding exhibited in the charge 
density. The Fe-d states dominate the structure at the 
Fermi level. In contrast the Nb-d character contributes 
far less structure to the overall density of states. This re- 
flects the role of Nb as an electron donor in the system. 

Clearly the Fermi level lies near a depression in the 
density of states which has been cited as important to 
the formation of a SDW at stoichiometry±H. The density 
of states at E F is 3.6 eV _1 FU _1 . In a rigid-band approx- 
imation the impact of doping a small amount of holes 
(e.g. Nb-doping) will raise the total density of states and 
promote the development of ferromagnetism as is seen in 
its doping phase diagram experimentally. However, if 
we take the same rigid-band approach to electron dop- 
ing, such as in the Fe-doped system, then the density 
of states is expected to fall. This is in conflict with the 
experimental phase diagram (see Fig. [I]), which shows 
the formation of ferromagnetism with small amounts of 
Fe doping. Therefore, effects beyond rigid-band shifts of 
Ep may be in operation and we attempt to explore these 



in the following section. 



III. DOPING DEPENDENCE OF THE 
FERROMAGNETIC INSTABILITY 

We may evaluate the proximity of the system to ferro- 
magnetism within Stoner theory^ 2 , where if N(Ep)I > 1 
then a material is unstable to ferromagnetism. I t has 
been showrP^l that if we assume a typical interactio n^ 3 * 24 * 
/ ~ 0.7 — 0.9 eV then we obtain a Stoner enhance- 
ment of (1 - NtEp )!)' 1 ~ 3 for stoichiometric NbFe 2 . 
Experiment ally^M however, comparing the measured 
susceptibility with the bandstructure value of N(Ep) in- 
dicates a Stoner enhancement (1 — N(Ef)I)~ 1 > 100. 
This suggests firstly, that the interaction scale may be 
unusually large in NbFe2, I > I eV. Secondly, very slight 
increases in the value of N(Ep) of order a few percent 
should drive NbFe 2 to a ferromagnetic instability. In this 
section we utilise doped supercell calculations to inves- 
tigate the evolution of N(Ep) with both hole (Nb) and 
electron (Fe) doping. We note that doped NbFe 2 sam- 
ples are disordered as is evidenced by the significant in- 
creases in the low temperature resistivities that occur 
upon doping. For these calculations we have employed 
the lattice parameters derived from experiment 7 and re- 
laxed all internal coordinates. We consider the non- 
magnetic density of states in order to apply the Stoner 
framework. 

First, we consider hole doping on the Nb-rich side of 
the phase diagram of Nt>i_ 2/ Fe2+ 2/ . We consider dop- 
ing corresponding to y = —0.065 Nb-rich Nb1.065Fe1.935. 
This state is formed by substituting one Nb atom for an 
Fe atom in a 2 x 2 supercell and lies beyond the for- 
mation of ferromagnetism in the experimental phase dia- 
gram. If we were to apply the rigid band approximation 
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to the density of states shown in Fig. |4j then we expect 
that for Nb1.065Fe1.935 the density of states at the Fermi 
level, N(E F ), will rise to 5.4 eV _1 FU _1 . Therefore the 
system would satisfy the Stoner criterion, N(E F )I > 1, 
for ferromagnetism. Clearly, since there are the two in- 
equivalent Fe(2a) and Fe(6ft) sites, then the Nb might 
be doped at either one. We have performed supercell 
calculations with the Nb atom substituted at both sites. 
Non-magnetic total energy calculation for the 2x2 su- 
percell indicated that for substitution at the Fe(6A) site 
the total energy is ~ 6.8eV/FU lower than that when we 
substitute at the Fe(2a) site. Therefore, the total energy 
strongly favours substitution at the Fe(6/i) site. 

In Fig.[5|a) we show the calculated density of states for 
Nb-rich Nb1.065Fe1.935, with the Nb substituted on the 
Fe(6A) site. Here we find that N(E F ) = 3.8 eV^FU" 1 . 
This is a smaller rise than that predicted by the rigid 
band approximation, but proabably still large enough to 
satisfy the criterion for the formation of ferromagnetism 
under Stoner theory. If we inspect the projected density 
of states for the Nb impurity and an Fe(6/i) site shown 
in the same figure, then we can see a mechanism for this 
result. The Fe(6/i) site contributes significant structure 
to the density of states near the Fermi level. As a result 
when we dope a Nb onto a Fe(6ft) site, not only are we 
removing valence electrons from the system, but we are 
also decreasing the contribution from Fe(6/i) sites to the 
density of states at the Fermi level. The importance of 
the structure of Fe-d states has been noted by Inoue et 
al in tight-binding calculations incorporating impurity 
atoms 25 . Consequently, the density of states at the Fermi 
level from our doped supercell calculation is lower than 
that expected in the rigid band approximation. 

Next we consider electron doping on the Fe-rich side 
of the phase diagram of Nbi_ 2y Fe2+ 2/ . We consider dop- 
ing corresponding to y = 0.065 Fe-rich Nbo.935Fe2.065- 
This state is formed by substituting one Fe atom for 
a Nb atom in a 2 x 2 supercell and is well within the 
ferromagnetic region of the experimental phase diagram. 
However, according to the rigid-band approximation, this 
doping would actually produce a small fall in the density 
of states at the Fermi level to N(E F ) = 3.5 eV _1 FU _1 . 
This is due to the fact that the Fermi level effectively 
moves higher in energy in Fig. [4] and as a result the 
density of states falls deeper into the trough near the 
Fermi level. Therefore, in the rigid band approximation 
the Stoner criterion would not be satisfied and ferromag- 
netism would not result. In Fig.[5|b) we show the density 
of states from our supercell calculation of Nbo.935Fe2.065- 
Critically, the density of states at the Fermi level actu- 
ally rises to N(E F ) = 4.5 eV^FlT 1 . This is due to the 
fact, that despite the addition of valence electrons, the Fe 
dopant contributes increased structure to the density of 
states at E F . To illustrate this we show in Fig. [5jb) the 
projected density of states for the Fe dopant at the Nb 
site. We see that, unlike Nb, this Fe dopant contributes 
a form to the density of states that is very similar to that 
of an Fe(6/i) site. 
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FIG. 5: (Color online) Non-magnetic calculated DOS for (a) 
Nb-rich Nb1.065Fe1.935 and for (b) Fe-rich Nbo.935Fe2.065- The 
insets show the same data on an expanded scale around the 
Fermi level. The Fermi level is set to eV. 

Also, we have made a simple approximation to the dis- 
order by applying a gaussian broadening of width 41meV 
to our density of states to reflect the decreased life- 
time of electrons in the states near the Fermi level due 
to disorder. This is calculated from the measur ed low 
temperature resistivity in a Drude approximatiorP^^. 
In Nb-rich Nb1.0e5Fe1.935 we find that N(E F ) = 4.6 
eV _1 FU _1 . The reason for the large increase compared 
to N(E F ) = 3.8 eV^FU" 1 found without this broaden- 
ing may be ascertained from the inset to Fig. [5|a). The 
Fermi level lies in a sharp dip in the density of states and 
therefore disorder induced broadening is likey to smear 
out this structure and increase N(E F ). This suggests 
that disorder may be important to the development of 
ferromagnetism upon Nb doping. In contrast N(E F ) 
is found to be insensitive to the broadening for Fe-rich 
Nbo.935Fe2.065- The influence of disorder may also ex- 
plain discrepancies between experimentally determined 
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phase diagrams where some worlP shows the develop- 
ment of ferromagnetism for lower Nb dopings than in 
more recent studied^. This indicative result invites fur- 
ther work on the description of disorder with approaches 
such as the Coherent Potential Approximation^HSI] A 
further consideration in such work may be how the inter- 
action, /, changes with disorder. 

These findings illustrate several important points 
about the evolution of the temperature-doping phase di- 
agram. 1) A rigid-band approximation is not sufficient 
to understand the evolution to ferromagnetism, particu- 
larly for electron doping. 2) The type of dopant will be 
critical to the character of this phase diagram. This is es- 
sentially because the fine structure of this phase diagram 
is strongly affected by the contributions to the density of 
states at the Fermi level from Fe-d states. 3) The pres- 
ence of disorder may be important to the development of 
ferromagnetism upon Nb doping. 

We have also undertaken calculations for dopings of 
y = ±0.0325 using 2x2x2 supercells. These calcu- 
lations deliver results with the same qualitative conclu- 
sions. For example we considered doping correspond- 
ing to y = 0.0325 Fe-rich Nbo.9675Fe2.0325- This state is 
formed by substituting one Fe atom for a Nb atom in a 
2x2x2 supercell. Applying the rigid-band approxima- 
tion would result in a fall from N(E F ) = 3.6 eV^FU" 1 
in stoichiometric NbFe 2 to N(E F ) = 3.0 eV _1 FU _1 . Our 
supercell calculation for Fe-rich Nbo.9675Fe2.0325 indicates 
N(Ep) = 4.1 eV _1 FU _1 and supports the formation of 
ferromagnetism as shown in the experimental phase dia- 
gram of Fig. [l] 



IV. FERMI SURFACE AND MAGNETIC 
INSTABILITIES 

To further our understanding of the complex evolu- 
tion of the magnetic phase diagram with doping we here 
consider the magnetic interactions active in stoichiomet- 
ric NbFe2, which has experimentally been interpreted as 
having SDW character. The previous work of Subedi et 
at. has considered the interaction of different magnetic 
orders within the unit celt^ under the local spin density 
approximation (LSDA). 

The Fermi surface of metallic materials is crucial in 
the understanding of their properties. In Fig. [6] we show 
the calculated non-magnetic Fermi surface for stoichio- 
metric NbFe2 which agrees well with previous worlP^. 
We reproduce this Fermi surface to emphasise the im- 
portance of features that may be important to the mag- 
netic instabilities of the system. Five sheets make up 
this Fermi surface which is strongly 3D. We consider the 
non-magnetic Fermi surface because we are interested in 
the magnetic instabilities as we approach the QCP where 
the static moment and spin splitting of the Fermi surface 
approach zero. 

We proceed by calculating Xo(q) using the technique 



of Romberg integration. 



Xo(q) = ^M, 

n, k 



2 

k,k+q 



/n,k «/n,k+q 
£-n,k+q £-n,k 



(i) 



Here, the sum is over all bands up to leV above the Fermi 
level, n, and wave vectors k. M are the matrix elements 
and / is the Fermi occupation function. No increased 
lifetime broadening was applied. In Fig. [6] (f) we show 
the susceptibility for q parallel to [010] and for q parallel 
to [001]. We have included only intraband contributions 
for which we may to a good approximation set the ma- 
trix element M = 1. Both susceptibilities rise steeply 
away from q = 0. Along [010], the susceptibility peaks 
at a small wave vector, q « [0 0.075 0]7r/b, and then 
falls monotonically until the zone boundary is reached. 
Along [001] the susceptibility rises to a higher level than 
that in-plane « 40% above that at q = 0. It then levels 
off to an essentially constant value until the zone bound- 
ary is reached. Enhancements at finite q are consistent 
with previous workH From these calculations it appears 
likely that finite q instabilities, particularly along [001], 
may be important in this system. 

We note also, that in the absence of spin orbit effects 
the matrix element for interband contributions is M = 
at q = due to orthogonality, then the inclusion of inter- 
band contributions will enhance the finite q instabilities 
shown in Fig. [6jf ) . 

Inspection of the contributions of each Fermi surface 
sheet shows that the rising value of Xo(q) along [001] 
is due largely to bands 82 and 83. This results from 
the presence of the very flat sections of surface that they 
possess. These flat sections nest strongly as they are 
translated by the wavevector q. This behavior is very 
similar to the large finite q instability found due to the 
flat sides of the cube-like Fermi surfaces in metamagnetic 
TiBe^. 

We also compare competing magnetic ground states by 
performing total energy calculations. Our total energy 
calculations consider several magnetic states including 
1) Non-magnetic, 2) Ferromagnetic, 3) AFM1 which is 
commensurate along the c-axis as suggested by the cal- 
culation of Xo(q)> 4) AFM2 which is q = antiferromag- 
netism of the type found in the C14 Laves material TiFe2 
and 5) AFM3 which is the lowest energy state found by 
Subedi et a/M We show the spin configurations associ- 
ated with some of these orders in Fig. [7J 

For the energy evaluation we have used 10,000 k-points 
in the Brillouin zone and the GGA-PBE exchange corre- 
lation functional. As shown in Table Ull we find that the 
lowest energy state corresponds to q = antiferromag- 
netism of the type found experimentally in the C14 Laves 
material TiFe2 33 . However, our results indicate that the 
Fe(2a) possesses a moment, while the experiments on 
TiFe2 indicate no ordered moment. In this calculation 
the second lowest energy state is then the ferrimagnetic 
state found to have the lowest energy by Subedi et a/P. 
The use of the LSDA correlation functional in that study 
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FIG. 6: (Color online) In (a)-(e) we show the non-magnetic Fermi surface of NbFe2. The flat sections of Fermi surface in 
bands 81 and 82 show potential for nesting. In (f) we show the non-interacting susceptibility, Xo(q), for NbFe2 along both the 
[010] and [001] directions. The large flat enhancement in Xo(q) along the [001] direction indicates the potential importance of 
large q magnetic instabilities. The flat sections of the Fermi surface of the bands 82 and 83 are the dominant contribution to 
the susceptibility at large q. A similar phenomenon is observed due to the cubic shaped Fermi surface in TiBeJ^. 
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FIG. 7: (Color online) Schematic diagram of the spin con- 
figuration for (a) AFM2 q = antiferromagnetism and (b) 
AFM3 the lowest energy state found by Subedi et alW^. Only 
Fe atoms are shown and the 2a and 6h rows are as indicated. 



results in a different configuration of spins for the lowest 
energy state. 

For the SDW order AFM1 corresponding to the peak 
in Xo(q) along the c-axis we find a lower energy than 
that obtained for ferromagnetism, but not lower than 
the configurations within the unit cell. The ground state 



TABLE II: Comparative energetics of magnetic states. For 
the ferromagnetic state the total magnetic moment is 2.34/is 
per formula unit. In the ferromagnetic state the Nb is found 
to also carry a small moment ~ 0.2/jb which is antiferromag- 
netically coupled to the Fe moments. 



Energy (meV/FU) 


Fe(2a) fi(fi B ) 


Fe(6/i) /jl(/jlb) 


NM 











FM 


-84 


1.1 


1.5 


AFM1 


-96 


1.6 


1.4 


AFM2 


-104 


1.6 


1.5 


AFM3 


-97 


1.8 


0.98 



calculation for these SDW configurations also produces 
large moments and suggests that the energy gain is again 
dominated not by the Fermi surface nesting, but by the 
Hund's coupling. 

Importantly, the ordered moment below the SDW 
transition (< 0.02/ie) represents only a small fraction 
of the fluctuating momentPl Magnetic order is then 
expected to have a minor effect on the magnetic fluc- 
tuation spectrum and hence on the low temperature 
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TABLE III: Results from spin fluctuation theor^ in 3D 
shown alongside the experimental results for NbFe2 at its 
quantum critical poinds Results are shown for the standard 
Fermi liquid (FL), a system near a finite q antiferromagnetic 
instability (AFM) and a system near a ferromagnetic insta- 
bility (FM). 



Property 


FL 


AFM 


FM 


NbFe2 


Specific Heat (C) 


T 


T 


Tln(T*/T) 


T\n(T*/T) 


Resistivity (p) 


ji2 






ji3/2 



properties, in particular the non- Fermi liquid low 
temperature exponents. This is evidenced by the small 
size of the heat capacity anomaly at Tjy ~ 10 K even in 
stoichiometric NbFe2- Therefore, the presence of finite 
q magnetic fluctuations may have an important role in 
producing the non-Fermi liquid properties. 

From spin fluctuation theory, it is known that at an 
antiferromagnetic (q ^ 0) quantum critical point the 
expected temperature dependence of the resistivity is 
p oc T 3 / 2 in the presence of sufficiently strong quenched 
disordei^, while the leading order term in the specific 
heat is linear in T. In contrast at a ferromagnetic crit- 
ical point (q = 0), p oc T 5 / 3 and the specific heat di- 
verges logarithmically. As shown in Table |III[ exper- 
iments on NbFe2 at its quantum critical point might 
suggest the presence of both ferromagnetic and antifer- 
romagnetic spin fluctuations if we were able to treat these 
modes separately. Alternatively, NbFe2 may be entering 
into a new regime such as that found in the high pres- 
sure region of the helimagnet MnSPn It is possible that 
the presence of q = and finite q instabilities in this 
material may combine to produce the unusual non- Fermi 
liquid behavior. The prominence of such instabilities in 
NbFe2 is confirmed by both our calculation of the sus- 
ceptibility in Fig. |6|f) and the total energy calculations 
in Table [TT1 

The interplay of Fermi surface nesting along with other 
magnetic couplings in this system is shared with the Fe- 
based superconductors where nesting is important along 
with other band structure effects. The competition of 



these effects has been shown by YildirinP^to produce the 
structural phase transition in the Fe-based superconduc- 
tors. In the NbFe2 system they may produce intriguing 
non-Fermi liquid properties. 



V. SUMMARY AND DISCUSSION 

The bonding structure of the Fe-cages in NbFe2 is 
likely to be instrumental in producing the magnetic in- 
terations that dominate the formation of its magnetic 
ground state. This bonding produces inequivalent char- 
acter in the kagome triangles in the Fe(6ft) layers. 

We demonstrate that non-rigid band effects must be 
considered in order to understand the evolution of the 
ground state as a function of doping. The contribution 
of Fe-d states to the density of states at the Fermi level is 
shown to be critical to the production of a Stoner instabil- 
ity upon electron doping with Fe impurities. Therefore, 
it is expected that the nature of the electron dopant will 
be crucial to the evolution of the phase diagram. Also, 
the presence of disorder is likely to be important to the 
formation of ferromagnetism upon Nb doping. 

We find that NbFe2 is close to both small and large q 
magnetic instabilities. The non-interacting susceptibility 
indicates the potential importance of q ^ instabilities. 
The electronic structure is highly three dimensional, but 
the possibilities for nesting of quasi-lD sections of the 
Fermi surface may provide crucial magnetic instabilities. 
Total energy calculations indicate that spin configura- 
tions within the unit cell produce the lowest energy state. 

The unusual non-Fermi liquid critical exponents of 
NbFe2 may arise from the combination of q ^ and 
q = instabilities. Further investigation as to how 
these instabilities may interplay will be necessary both 
theoretically and experimentally. Furthermore, the 
presence of such a large phase space for fluctuations is 
likely to complicate the formation of superconducting 
Cooper pairs in either the spin singlet or triplet channels. 

The authors would like to acknowledge useful discus- 
sions with J.B. Staunton. 
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